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ABSTRACT 

We review the accuracy of existing Fe XVII X-ray line emission models by 
comparing them with an extensive analysis of Chandra high energy transmission 
grating (HETG) observations of stellar coronae. We find significant discrepan- 
cies between most theoretical predictions and observations for at least some of 
the intensity ratios involving the six principal Fe XVII lines, 3C (15.01 A), 3D 
(15.26 A), 3E (15.45 A), 3F (16.78 A), 3G (17.05 A), and M2 (17.10 A). We sug- 
gest that the main problem of most previous theoretical studies to their inability 
to fully include electron correlation effects in the atomic structure calculations, 
while any deficiencies in the scattering approximation methods are of minor im- 
portance, regardless of it being close-coupling (CC) or distorted- wave (DW). An 
approximate method based on the many-body perturbation theory and DW ap- 
proximation is proposed to include such correlation effects in the calculation of 
collisional excitation cross sections. The results are shown to agree with coronal 
observations and laboratory measurements better than most previous theories. 
Using the new atomic data, we then investigate the electron density sensitivity 
of the M2/3G intensity ratio and provide an improved density diagnostic tool for 
astrophysical observations. 

Subject headings: atomic data — atomic processes — line: formation — X-rays: 
general 



1. Introduction 



There are several widely known problems in the theoretical modeling of Fe XVII X- 
ray emission lines, making it difficult to use them as reliable diagnostic tools for electron 
density, temperature, plasma opacity, and iron abundances in X-ray astronomy. The first 
problem is related to the intensity ratio of the 2p 3 / 2 -^d5/2 transition at 15.26 A, gener- 
ally referred to as 3D, to the 2pi/2-3d 3 /2 transition, or 3C. Most th e oretical predictions of 



this ratio are smaller than astrophysical observations (IXu et al.ll2002l ; iBrinkman et al.l 12000 
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farown et al.lll99 j 



Canizares et all l2000t iMewe et al.ll200ll; iBehar et al.l l200ll ) and laboratory measurements 



2001bl ; iBeiersdorfer et al.ll200ll ). This discrepancy has been used in the 



past to infer plasma opacity as the stronger 3C line is more susceptible to resonant scat- 
tering (IWaljeski et al.lll994j ; ISaba et al.lll999l ). However, this diagnostics depends critically 
on knowing the optically thin limit of the 3D/3C ratio. The second problem concerns the 
intensity ratio of the 2p-3s (hereafter 3s) transitions, including 3F (16.78 A), 3G (17.05 A), 
and M2 (17.10 A), to the 2p-3d (hereafter 3d) transitions, including 3C, 3D, and a much 
weaker 3E (15.45 A). The astrophysical observa tions and laboratory measurements are also 



signi ficantly larger than most theoretical values (jPhillips et al.lll999l ; IBeiersdorfer et al.ll2002 



2004 ) for this ratio. Finally, the theoretical predictions of the intensity ratio of the two 2p-3s 
transitions, M2 to 3G, have also been shown to be significa ntly smaller than astrophysical 
observations of low density plasmas (INess fe Schmittl 120051 ). The M2/3G ratio is density 
sensitive, and have b een used to infer electron densities in magnetic cataclysmic variables 
(IMauche et al.ll200ll ). However, the inability of reproducing the low density limit of this 
ratio in theory introduces additional uncertainties in the density determination. 

There have been claims that these modeling problems have been solve d with large scale 
close- coupling calculat ions of the electron collisional excitation data of (IChen fe Pradhan 
20021 ; IChen et al.l 120031 ). The new theoretical values of certain line ratios have been shown 
to agree with laboratory measurements and a few astrophysical observations to wit hin 10%. 
However, the cross section measurements of the 3C and 3D lines ( Brown et al 2006h indicate 
that such apparen t agreements are rather fortuitous. The more recent work of IChenl (120071 ) 
and IChenl (l2008bl ) revised their earlier work by including pseudo states and the 3d 2 pair ex- 
citation configuration in the target wave function expansion and obtained 15-20% reductions 
in t he 3C cross s ections compared with IChen et al.l (120031 ). while those of 3D changed lit- 
tle. IChenl (j2008al ) further suggested that the theoretical radi ative recombination (RR) cross 
sections into the M-shell of Fe XVI used for normalization in iBrown et al.l (120061 ) are under- 
estimated by 25%, leading to smaller measured 3C and 3D cross sections. In this paper, we 
provide a thorough assessment of the accuracies of existing atomic calculations on Fe XVII 
by comparing them with a large sample of Chandra observations of stellar coronae obtained 
with the high energy transmission grating (HETG) spectrometer. We conclude that most 
of the existing predictions do not offer satis factory expla nation of the observed line ratios. 
The recent calculations of IChenl (120071 ) and I Chen! (l2008bl ) provide the best agreement with 
the laboratory 3C and 3D cross sections. However, a detailed comparison of these data with 
the Chandra observations cannot be done, as the full results of the calculations including 
al l levels neces sary for col li sional radiative modeling are not available. The major drawback 
of IChenl (120071 ) and I Chen! (j2008bl ) is that it is not clear that the accounting the target cor- 
relation effects has fully converged. In these calculations, the largest improvement comes 
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with the inclusion of 3d 2 pair excitation configuration. However, it is well know that such 
pair correlations are notoriously slowly convergent. In this paper, we propose a many-body 
perturbation correction procedure for the excitation cross section calculations, taking all 
singly and doubly excited correlations into account, but treating their eff ects o n excit ation 
cross section s approximately. We show that the results agree well with IChenl ( 120071 ) and 



Chen! (j2008al ) for the 3C and 3D cross sections, and all Fe XVII line ratios are brought into 



agreement with the Chandra observations. 



2. Review of existing theoretical calculations 



Th eoretical study of Fe XVII X-ray lines and their applications in astronomy have a long 
history. iBely fe Belyl (119671 ) obtained t he excitation cross sections of 2p 5 3/ con figurations us- 
ing the Coulomb-Born approximation. lLoulergue fe Nussbaumerl (119731 . 11975T ) developed col- 
lisional radiative models of Fe XVII and Ni XIX with DW collision stre n gths, and compared 
the re sultin g line intensities with s olar observations. ISmith et al.l (119851 ). iRaymond fe Smith 
(119861 ). and lGoldstein et al.l (Il989l ) pointed out the importance of resonan t excitation and its 



role in determining the temperature dependence of Fe XVII line ratios. iBhatia fe Doschek 



tion. 



(11992) presented a comple te set of Fe XVII line intensities obtained with the DW approxima- 



Cornille et al.l (119941) studied the Fe XVI satellite transitions that blend with Fe XVII 



lines. iMohan et aLT (jl997t ) investigated the effects of relativistic term coupling effects of the 



Fe XVII collision strengths using the Breit-Pauli i?-matrix approximation. 



Beiersdorfer fe Wargelinl (119941 ) built a flat crystal spectrometer on the electron beam 



ion trap (EBIT) of the Lawrence Livermore National Laborato ry (LLNL), and co nducted 
the first laboratory measurement of Fe XVII X-ray line emission. iBrown et al.l (119981 ) carried 
out a systematic investigation of the Fe XVII lines with the same spectrometer and EBIT 
facil ity, and showed th a t most prey ious theories significantly underestimate the 3D/3C ra- 



tio 



Beiersdorfer et al.l (120021 . |2004j) showed that the 3s /3d ratio is also underestimated in 
theory by large factors as compared with both EBIT and Tokamak measurements. Inde- 



pendent laboratory measurements by lLaming et al.l (120001 ) with the National Institute of 
Technology and Standards (NIST) EBIT obtained 3D/3C ratios that are consistent with 
those from the LLNL EBIT, but reported 3s/3<i ratios very different from those measured 
at the LLNL EBIT and concluded that the 3s /3d ratios are consistent with modern theoret- 
ical calculation s. These laboratory work prompted new theoretical stud ies to resolve these 
discre pancies. IBhatia fe Doschek! (120031 ) improved the earlier work of IBhatia fe Doschek 



(119921 ) by including mo r e con figurations and computing collision strengths in a wider energy 



range. 



Doron fe Beharl (120021 ) investigated the effects of resonant excitation, recombination 
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of Fe XVIII and inner-shell ionization of Fe XVI as line formation processes using the Hebrew 
University Lawrence Livermore Atomic Code (HULLAC), and concluded that such a multi- 
ion approach increases the 3s /3d ratios, bringin g th e m in b etter agreements with laboratory 



measurements and astrophysical observations. iGul (120031 ) reached simila r conclusions and 
also showed that recombin atio n processes are important for other iron ions. I Chen fe Pradhan 
( 120021 ) , I Chen et al.l (120031 ) and I Chen fe Pradhanl (120051 ) carried out a large scale Breit-Pauli 
.R-matrix calculation for electron collisional excitation of Fe XVII, and claimed that their 
results resolve the discrepancies between theories, laboratory measurements, and astrophysi- 
cal ob servations. An in dependent close-coupling calculati on with the Dirac a tomic .R-matrix 
code (lLoch et al.l 120061 ) confirms some of the findings in IChen et al.l (120031). bu t foun d re- 
maining discr epancies between theories and measurements. Most recently, IChenl (120071 ) and 
Chenl ( I2008bl ) revised their earlier work by including pseudo states and the 3d 2 pair excita- 
tion configuration i n the target wavefunciton expansion, and obtained significantly lower 3C 
cross sections than lChen et al.l (120031 ) . 



In the present paper, we concentrate on these newer theoretical work. A comparison 
of the Fe XVII line ratios from these calculations are shown in Figure H} alo ng with the 



predic tions from the widely used astrophysical plasma emission code (APEC) of ISmith et al. 



(120011 ). which is based on DW collisional data without resonances. All line ratios involved 
are calculated at the low electron density limit of < 10 10 cm" 3 . It is immed i ately clear from 
this comparison that t he predictions from A P EC a nd iBhatia fe Doscheld (12003!) are very 



different from those of iDoron fe Beharl (120021 ). IGul (120031 ). and lLoch et all (120061 ). These 



differences can be attributed to the lack of resonant excitation contributions in APEC and 



Bhatia fe Doschek! (120031 ). On the other hand, the results of IDoron fe Beharl (120021 ) and 
Gul ( 120031 ) are relatively close to each other both in the temperature dependence of the 
line ratios and their magnitudes, which is not surprising since the two calculations employ 
essentially the same physics and appro xima tion s. The relatively large difference in the 3E/3C 



ratio between IDoron fe Beharl (120021) and 



treatment o: 



Gul (120031) can be attributed to the incom plete 



reso nant excitation in IDoron fe Beharl (120021 ) . T he results o f IGul (120031 ) and 



Loch et al.l (120061 ) are also in very good agreement, although lLoch et al.l (120061 ) generally 



produce a steeper tempera ture depend e nce o f the line ratios, especially for the 3s/3<i ratios. 
This is due to the fact that lLoch et al.l (120061 ) do not include recombination contributions to 
the line intensities. Recombination preferentially enhances the 3s transitions and has larger 
contributions at higher temperatures due to the higher fractional abundance of Fe XVIII 
relative to Fe XVII. One sur prising fact i n this comparison is the relatively large differ ences 



between the 3F/3C ratios o flChen et a 



the corresponding ratios of lLoch et al 



J2003I) . 3s/3d ratios of lChen fe Pradhanl J2OO5I) . and 
(I2006I ). In fact, the ratios of IChen et al.l J2OO3I ) and 



Chen fe Pradhanl (120051 ). which include resonant excitation, are similar to the DW values 
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without resonant excitation. lLoch et al.l (120061 ) showed that the differences in the 3C and 
3D cross sections in the two i?-matrix calculations are small. Therefore the large differences 
we note in this comparis on are due to the differences in the 3s intensities. The updated 3C 
and 3D cross sections of Chen J2007T ) give larger 3D/3C ratio s than all other calculations. 
Because the full results of the new calculations of IChenl (120071) are not availa ble, we scaled 
the 3F/3C and 3s/3d ratios of Ichen et all J200J and lchen & Pradhanl t)2005h to reflect the 
updated 3C and 3D cross sections, b ut k e eping 3s in tensities unchange d. These scaled ratios 
for 3F/3C become close to those of iGul (120031 ) and lLoch et al.l (120061 ). but those for 3s/3<i 
are still significantly smaller. 



3. Chandra stellar observations of Fe XVII spectra 



In order to assess the accuracies of various calculations, we compare them with an exten- 
sive set of stellar coronal observations with the HETG spectrometer aboard Chandra X-ray 
observatory. Since its launch, Chandra has observed a large sample of stellar coronal sources 
with the high resolution HETG spectrometer. Table [1] lists the sources in our selection. 
The data were retrieved from the Chandra archive, and reprocessed with CIAO version 3.4 
and CALDB version 3.4. Some sources are observed multiple times, and we combine the 
individual observations to improve the statistics. There are 11 sources that do not provide 
decent counting statistics individually, and we have also combined them to produce a single 
composite observation referred to as "Stack" in Table [TJ Of the 24 targets (after grouping of 
multiple observations and weak sources) in our sample, Capella provides the best statistical 
quality because it is an calibration target, and observed many times. 



We follow the differential emission measure (DEM) analysis method of iGu et al.l (120061 ) 
to characterize the coron al propert y of these sources. Using the APEC database modified 
by the line emissivities of IGul ( 120031 ) and wavelengths of IGul (120051 ). we reconstruct the DEM 
distribution and elemental abundances of the coronal plasma by jointly fitting the ±1 orders 
of medium and high energy grating (MEG and HEG) spectra. The resulting DEM and 
comparison of theoretical and observed spectra in the 12-18 A region is shown in Figures 
[2] and [3] for the composite observati o n of C apella. The DEM is peaked near 10 6 9 K, which 
is similar to the result of IGu et al.l (120061 ) using a single observation. To characterize the 
temperature of the Fe XVII emission regions, it is more appropriate to examine the DEM 
weighted by the Fe XVII emissivities, which is also shown in Figure [21 It indicates that 
the Fe XVII emissivity weighted DEM has a narrower profile. In fact, the characteristic 
temperatures of the Fe XVII emission regions are always in the range of 10 6 ' 6 -10 6,9 K for all 
sources, even though some of them have significant high temperature tails beyond 10 7 K. 
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This is due to the fact that Fe XVII fractional abundance peaks near 10 6 6 K in collisional 
ionization equilibrium. 

The spectral model shown in Figure [3] already provides a reasonable fit to the observed 
spectrum. However, close inspection of the Fe XVII lines reveal significant discrepancies 
between them. Using the reconstructed spectrum as a basis, we obtain a bet ter fit to the 



Fe XVI I emi ssion by varying 3C, 3D, 3E , 3G, 3F, and M2 line intensities. iBrown et al. 
(j2001al ) and iBrickhouse fe Schmela (120061 ) demonstrated that the 3D line blends perfectly 
with an Fe XVI satellite transition, and another nearby Fe XVI transition at 15.21 A can be 
used to remove the Fe XVI contribution to the 3D line. Many of the sources in our sample 
indicate the presence of the 15.21 A line, providing a measure of contamination to the 3D 
line. The 15.21 A line intensity is extracted from the measured spectra and multiplied by a 
factor of 0.7, which is then subtracted from the apparent 3D intensity. The branching ratio 
of 0.7 used here is determined w ith the Flexible Atomic Code developed by the author, and is 



betwee n the values 0.83 used by IBrickhouse fe Schmelzl (120061 ) and 0.51 used by lBrown et al. 



(j2001al ). However, because the Fe XVI contamination of the 3D line is relatively small, the 
correction decreases the measured 3D/3C ratio only by ~ 10%. The differences in the 
branching ratios leads to 2-3 % differences in the 3D/3C ratio, which are smaller than the 
statistical errors for most sources. The derived line ratios for the 24 sources are listed in 
Table [2 ] The Fe XVII em issivity weighted average temperature, T e m, of the coronae are also 
shown. iNess et al.l (120031 ) reported the ratios involving the 3C, 3D and 3F lines for a large 
sample of coronal sources observed by C handra and XMM -Newton. The results obtained 
here are generally consistent with those of INess et al.l (120031 ) for the set of common sources. 



In comparing the measured and theoretical line ratios, the calculations should ideally 
be performed for the appropriate DEM distributions. In practice, we find that the ratios 
calculated with the DEM is very close to the ones calculated at a single temperature of T em . 
This is demonstrated in Figure H] for the 8 line ratios we investigate. In the rest of the paper, 
we therefore use the isothermal line ratios calculated at Tem for comparison. 

Figure [5] shows the comparison of measur ed and theoretical line ratios in detail. Be - 
cause the pre dictions of iDoron fe Beharl (120021 ) are generally similar to those of iGul (120031 ). 
and those of iBhatia fe Doschekl (120031 ) are similar to those of APEC, we omit those from 
Figure [5] to reduce clutter. It is clear that none of the theories can explain all line ratios 
satisfactorily. All theoret ical ratios are below the measured ones, except for 3D/3C, where 

agree with the measurem ents, 3E/3C and 3 F/3C, where the values 



the results of Chen 



of lLoch et all ( 200 



Gul (120031) . and the results of IChen et al.l (120031 ) scaled for the new 3C 



Chen J2007h agree with the measurements, and M2/3G, where the values of 



cross sections of 

APEC agree with the measurements. However, it should be pointed out that the measured 
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3E/3C ratio have large statistical uncertainties. The APEC ratios of M2/3C and 3G/3C are 
very different from the measurements, and the agreement for M2/3G is therefore fortuitous. 
Despite the claim of excellent agreements bet ween their theoret i cal line ratios, laborator y 
mea surements and ast r ophys ical observations in I Chen &: Pradhanl (120021 ). IChen et al.l (120031 ). 
and IChen fe Pradhanl (120051 ) , we find their theoretical values are significantly smaller than 
the Ch andra measuremen ts for the 3F/3 C and 3s /3d ratios, even more so than the calcula- 
ti qns of Loch et al.l (120061 ) and iGul (120031 ) . After scaling for the new lower 3C cross sections 
of IChenl tooll the 3 s /3d r atios are still smaller than the Chandra observations. Since the 



Chen et al.l (120031 ) and subsequent work are not publicly available, we are not able 



data from 

to calculate other line ratios for thermal plasmas for a more detailed comparison. 



4. Distorted-wave with many-body perturbation correction 



The fact that the calculations of IGul (120031 ) and lLoch et al.l (120061 ) agree with each 
other very well except for the differe nces i n the recombination co ntributions i ndicat es that 
the two theoretical methods, DW in IGul (120031 ) and i?-matrix in lLoch et al.l (120061 ) are of 
comparable quality for Fe XVII. The differences between measured and calculated line ratios 
are not due to deficien cies in the sca t tering approximations. The cross section measurements 
for 3C and 3D lines of iBrown et al.l (120061 ) demonstrated that theories tend to overestimate 
the 3C cross sections more than 3D. The general trend in the comparison of theories and 
Chandra measurements indeed support t his picture . If t he theoretical 3C cross sections were 
to decrease by ~ 20%, the predictions of IGul (120031 ) and lLoch et al.l (120061 ) would agree with 
the observations much better. 3C is a dipole allowed transition, its collisional excitation cross 
section is nearly proportional to the oscillator strength, which is determined solely by atomic 
structure calculations. Therefore we suspect that the root problem in the previous theories 
of Fe XVII line intensities is the incomplete treatment of electron correlation effects in the 
target ion. This conclusion is supported by t he comparison of various oscillator strengths 
calculations for 3C and 3D listed in Table 2 of lChen et al.l (120031). Of these ca lculations, the 
many-body perturbation theory (MBPT) treatment of ISafronova et al.l (120011 ) is likely to be 
the most accurate in describing electron correlation effects in Fe XVII, which a lso yie l ds the 
sma llest oscillato r strength for the 3C transition. The recent calculations of IChenl (120071 ) 
and IChenl (l2008bl ) also showed that by improving the target wavefunction expansion, the 3C 
cr oss sections are r educed by 15- 20%. However, it is not clear that the target description 



m 



Chenl ( 120071 ) and IChenl (l2008bl ) has fully converged, given that the largest change comes 



with the inclusion of the 3d 2 pair excited configuration, and such pair correlations are know 
to be notoriously slowly convergent. 
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We have developed an independent MBPT method within the framework of FAC for 
general open-shell ions, a nd have u sed it in calculating the wavelengths of iron and nickel 
L-shell X-ray transitions jGull2005h . We recently extended the method to also calculate 



ra- 



diative transition rates and oscillator strengths. The 3C and 3D weighted oscillator strengths 
calculated in our MBPT method are 2.17 and 0.62, respectively. This compares to 2.50 and 
0.62 calculated by the DW approximation. These independent MBPT calculations therefore 
suggest that the 3C cross sections ar e likely to be overestimated by as much as ~ 10-20 % 



in previous non-MBPT calculations. IChenl (120071 ) gave the weighted oscillator strengths of 



3C and 3D to be 2.25 and 0.635, respectively, which agree with our MBPT results very well. 
This indicates the near convergence of their accounting of Fe XVII target correlation effects. 

However, the extension of the MBPT method to treat collisional excitation processes is 
rather difficult beyond the plane-wave Born approximation. Here we propose an approximate 
procedure for including; MBPT corrections in the calculation of DW excitation cross sections. 



Bar-Shalom et al.l (119881 ) showed that the DW collision strengths can be factorized to contain 



angular and radial integrals as follows 



k a o a i 



where a, (3 denote bound orbitals making the transition, Z k are angular factors, and Q k are 
radial integrals involving partial-wave sum 

Q k (a a 1 ;f3 /3 1 ) = X)^~ lp *( K ° /Ci; a ^i)P k {^K X ; fopi), (2) 

and 

P k (K QKl ; « «i) = X Wo; «i«i) + y2(-l) k+t [k] l 3ao 31 ] ) XWo; Riai), (3) 

t I Jo 3 ax k J 

where k represents a continuum orbital and X k are the usual Slater integrals. The first term 
in P k is the direct contribution, and the second term is the exchange contribution. In the 
Coulomb-Bethe approximation, the direct term can be written as 

P£b(« /«i; a ax) = M fc (a ai)-Rfc(«o«i), (4) 

where 

M fc (a ai) = <a \\C k \\ai> J(P ao P ai +Q ao Q ai )r k dr 
Rk(^oHi) = [ {Pk q Pkx + Q K0 QKx)zk+T dr ' ( 5 ) 
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where C k is the normalized spherical harmonic, P a and Q a are the large and small com- 
ponents of the Dirac radial wavefunctions of orbital a. Note that depends only on the 
target orbitals and is proportional to the electric multipole integral of rank k. This integral 
is the same as the one used in radiative transition rate calculations. Therefore, the ratio 
of the MBPT corrected to the conventional DW direct excitation contributions is the same 
as the ratio of the MBPT corrected to uncorrected oscillator strenghts. We calculate such 
ratios from the radiative transition rates with and without MBPT effects, and use them as 
scaling factors to correct the direct excitation contributions to P k , and obtain the MBPT 
corrected collision strengths in the DW approximation. 



Following the procedure outlined above, and the three- ion model of iGul (120031 ). we 
recalculate the Fe XVI I line ratios and compare them with the Chandra observations in 
Figure EJ The results of lLoch et al.l (120061 ) are shown again for easy comparison. It is obvious 
that our MBPT corrections bring the calculated line ratios into much better agreement with 
Chandra observations. The M2/3C and M2/3G ratios seem to be still slightly below the 
observed ones, but show significant improvements over the existing theories. 

Finally, we compare the MBPT corrected effe ctive formation cros s sections of the 3C 
and 3D transitions with laboratory measurements of iBrown et al.l (120061 ) in Figured Our ef- 
fective cross sections include contributions of resonant excitation and radiative cascades from 



highe r level s. The cross sections of IGul ( 120031 ). I Chen et al.l (120031 ). lLoch et al.l (120061 ). and 
Chen! ( 120071 ) are also shown for comparison. The i?-matrix calculations include resonances, 
but do not taken int o ac c ount radiative cascades. Clearly, the MBPT corrections leave the 
3D cross sections of IGul (120031 ) al most uncha n ged, but d ecrease the 3C cr oss sections by 
about 18%. The cross sections of IChen et al.l (120031 ) and lLoch et al.l (|20061 ) are all higher 
than the present MBPT corrected values for 3C. The cross sections of I Chen! (120071 ) are only 
slightly lower than the present calculations by ~5%, which likely represents the differences 
between different scattering approximations. However, the present MBPT corrected cross 
sections for both 3C and 3D are still about 15% higher than the laboratory measurements 
overall, while the 3D/3C ratios agree with each other very well. The energy dependence of 
the measured cross sections shown in Figure [7] were obtained by sweeping the electron beam 
energies and recording the energy dependence of the 3C and 3D line intensities. They are 
normalized to the RR emission into the M-shell of Fe XVI at a single energy point of 964 eV. 
There is significant statistical fluctuation in the measured data, but the overall discrepancies 
between the measured values and the present calculations are less than 15%, especially at the 
near thresho ld region bel ow 1200 eV. The discrepancies between the measured values and 
the results of lChenl (120071) are even smaller, at 10% level. If the s uggestion tha t the RR cross 
sections used in lBrown et al.l (120061 ) are underestimated by 25% (IChenll2008al ) i s correct, the 



new theoretical cross sections of 3C and 3D of the present paper and those of I Chen! (120071 ) 
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and I Chen! (l2008bl ) would have become smaller than the renormalized measurements. 

Since our new Fe XVII line ratios agree with Chandra observations, the lower experi- 
mental 3C and 3D cross sections would lead to lower cross sections for all other lines as well, 
and to inaccurate iron abundance d e termi nations if these lines are used as main constraints. 
The cross sections of iBrown et al.l (120061 ) were determined by normalizing the collisional 
excitation lines to the RR emission onto the M-shell of Fe XVI. The overall normalization 
uncertainties were estimated to be on the order of 10%. The differences between our new cal- 
culations and the measurements are therefore not significantly larger than the experimental 
uncertainties, especially if one takes into account that the more elaboate R-Matrix scattering 
method seems to further lower the cross sections by ~5%. 



5. Density diagnostics with the M2 line 



Due to the relatively slow radiative decay of the M2 transition, the M2 line is collisionally 
quenched at electron densities above 10 13 cm -3 . T he M2/3G ratio has been used to infer 
electron densities in magnetic c ataclysmic variables (IMauche et al.ll200ll ) and stellar coronal 
sources (INess fe Schmittl 120051 ). However, £J3] shows that the M2/3G ratios of all coronal 
sources with densities known to be below 10 12 cm -3 are significantly larger than previous 
theoretical pred i ctions , which is identi fied as a major problem in applying this diagnostic in 
Ness fe Schmittl (120051 ) . We note that INess fe Schmittl (120051 ) seem to indicate that the low 
density limit of M2/3G calculated with lGul f|2003n " data agree with the co ronal obse rvations. 
We have since confirmed that it is du e to some misinterpretation of the iGul (120031 ) data in 
calculating the M2/3G ratio, and the IGul (120031 ) ratio of M2/3G is in fact much lower than 
the observations. The present MBPT corrected M2/3G ratio, although appears to be still 
slightly lower overall, is at least marginally consistent with low density coronal observations. 

In Figure [U we show the density dependence of the emissivities of the six Fe XVI I 
lines at a temperature of 10 6,7 K in compa r ison w ith the predictions of lLoch et al.l (120061 ). 
The publicly released data from lLoch et al.l (120061 ) do not include forbidden transition rates. 
In order to derive the density dependence of line emissivities, we augmented their dataset 
with the forbidden transition rates from th e present calculatio ns. Two major differences are 
seen between the present work and that of lLoch et al.l (120061 ). namely, the smaller 3C (due 
to MBPT correction) and larger M2 (due to recombination contributions) in the present 
calculation, which helps to explain the observed coronal values. The densi ty dependence 
of M2/3G ratios in the present calculation and those of lLoch et al.l (120061 ) are shown in 
Figure [9] f or two electron tem peratures, 10 67 and 10 72 K. The differences between the present 
work and lLoch et al.l (120061 ) are even more pronounced at higher temperatures due to the 
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increasing importance of recombination contributions to the M2 line intensity. 



conclusions 



We have conducted a systematic study of the Fe XVII emission lines of stellar coronal 
sources observed by the HETG spectrometer aboard Chandra. At least some of the line 
ratios involving the six prin cipal Fe XVII lines ar e shown todisagree with most existing 



theories. The calculations of iDoron fe Beharl (120021 ) and iGul (120031 ). which include resonant 



excitation and recombination contributions in DW approximation pr ovide as good agre ement 



with data as the more sophisticated close-coupling calculation of lLoch et al.l (120061 ) with 



Chen & Pradhan 


(2002); 


Chen et al. 


(2003 





agreements with the Chandra observations for the 3F/3C and 3s/3<i line ratios. The source 
of discrepancies in t he two -R-ma trix calculations are not clear. The most recent Breit-Pauli 



i?-matrix results of iChenl (120071 ) gave the best agreement with laboratory measurements of 
3C and 3D cross sections, but a systematic comparison with Chandra line ratios cannot be 
carried out, as the full results of these calculations necessary for collisional ra diative modeling 
are not avail able. If w e assume that the 3s lin e intensities of IChenl ( 120071 ) are the same as 



Chen et al.l (120031 ) and I Chen fe Pradhan! (120051 ). and only correct the ratios for the updated 



3C and 3D cross sections, the resulting 3D/3C and 3F/3C ratios are in agreement with the 
Chandra observations, but the 3s/3<i ratios are still significantly lower than observed values. 

Based on the MBPT calculations of oscillator strengths, we suggest that the incomplete 
treatment of electron correlation effects of the Fe XVII target is the cause of overestimating 
the 3C cross section s in most p r evious theories, as is also suggested by the recent calcula- 
tions of IChenl (120071 ) and IChenl (l2008bl ). We propose a simple, MBPT based procedure to 
incorporate such correlation effects in the calculation of electron collisional excitation cross 
sections under DW approximation. The resulting Fe XVII line ratios are shown to agree with 
Chandra observations much better than previous predictions. The 3C and 3 D cross sections 
of the present work are only ~5% larger than the calculations of I Chen! (120071). Compar i son o f 
our MBPT corrected 3C and 3D cross sections and the measurements of iBrown et al.l (120061 ) 
suggests remaining discrepancies on the order of 15%, although systematic uncertainties of 
about 10% may explain part of this difference, and the simplied scattering approximation 
used here may be responsible for the rem aining discrepancies . We have shown that if the 
RR cross section s used for nor malization in lBrown et al.l (120061 ) are underestimated by 25%, 
as sugge sted by IChenl (2008a), t he new theoretical cross sections of the present work and 



those of lChenl (hoOTt ) and IChenl (l2008bl ) would have become smaller than the renormalized 
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measurements. 

We have presented the density dependence of the Fe XVII emissivities and M2/3G 
ratio, and compared them with previous calculations. The improved low density limit of 
M2/3G ratio provides greater confidence in applying this density diagnostics in astrophysical 
observations. 
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and on behalf of the National Aeronautics Space Administration under contract NAS8-03060. 
This work was performed under the auspices of the U.S. Department of Energy by Lawrence 
Livermore National Laboratory under Contract DE-AC52-07NA27344. 
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Table 1. Chandra HETG observations selected in the present sample. Targets with 
multiple exposures are combined to form a single observation. 11 weak sources are also 
grouped to form a composite observation for target 15, or "Stack" . 



Index 


Target 


Obsld 


Exposure 


1 


£ Uma 


1894 


70.9 


2 


AD Leo 


2570 


45.2 


3 


T SCO 


638 


59.2 


4 


EV Lac 


1885 


100.0 


5 


CC Eri 


4513 


89.5 


6 


HD 45348 


636 


94.6 


7 


AU Mic 


17 


58.8 


8 


44 Boo 


14 


59.1 


9 


VW Cep 


3766 


116.6 


10 a 


Capella 




297.9 




Capella 


57 


28.8 




Capella 


1099 


14.6 




Capella 


1100 


14.6 




Capella 


1101 


14.6 




Capella 


1103 


40.5 




Capella 


1235 


14.6 




Capella 


1236 


14.6 




Capella 


1237 


14.6 




Capella 


1318 


26.7 




Capella 


2583 


27.6 




Capella 


3674 


28.7 




Capella 


5955 


28.7 




Capella 


6471 


29.6 


11 


HD 93497 


3410 


57.0 


12 


AB Dor 


16 


52.3 


13 


TZ Crb 


15 


83.7 


14 


V824 Ara 


2538 


94.2 


15 b 


Stack 




711.7 




TY Pyx 


601 


49.1 




FK Com 


614 


41.4 




HD 68273 


629 


64.9 




II Peg 


1451 


42.7 




HD 206267 


1888 


34.1 




HD 206267 


1889 


39.5 




HD 111812 


1891 


130.2 




HD 223460 


1892 


95.7 




PROXIMA Cen 


2388 


42.4 




HIP 92680 


3729 


73.9 




V4046 Sgr 


5423 


97.7 


16 


ER Vul 


1887 


112.0 


17 


HR 1099 


62538 


94.7 


18 


P Ceti 


974 


86.1 


19 


Algol 


604 


51.7 


20 c 


AR Lac 




64.3 




AR Lac 


6 


32.1 




AR Lac 


9 


32.2 



(ks) 
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Table 1 — Continued 



Index 


Target 


Obsld 


Exposure (ks) 


21 


UX Ari 


605 


48.5 


22 d 


IM Peg 




192.1 




IM Peg 


2527 


24.6 




IM Peg 


2528 


24.8 




IM Peg 


2529 


24.8 




IM Peg 


2530 


23.9 




IM Peg 


2531 


23.9 




IM Peg 


2532 


22.5 




IM Peg 


2533 


23.7 




IM Peg 


2534 


23.9 


23° 


a Gem 




120.8 




a Gem 


5422 


62.8 




a Gem 


6282 


57.9 


24 


A And 


609 


81.9 



a Sum of the following 13 observations of Capclla. 

b Sum of the following 11 observations of different targets. 

c Sum of the following 2 observations of AR Lac. 

d Sum of the following 8 observations of IM Peg. 

°Sum of the following 2 observations of a Gem. 



Table 2. Fe XVII line ratios measured with the Chandra HETG observations. T em is the 
Fe XVII emissivity weighted average temperature of sources in units of 10 6 K. The numbers 
in parentheses are statistical uncertainties. 3s refers to the sum of 3F, 3G, and M2, and 3d 

refers to the sum of 3C, 3D, and 3E. 



Index Target T em (MK) 3D/3C 3E/3C 3F/3C 3G/3C M2/3C M2/3G 3s/3C 3s/3d 



1 


£ Uma 


5 


13 





423 





024) 





048 


;o 


009) 





780 





041) 





974 





048) 





969 


;o 


048) 





997 


;o 


063) 


2.723 





098) 


1.852(0 


062) 


2 


AD Leo 


5 


30 





360 





048) 





059 





018) 





633 





070) 





816 





089) 





814 




082) 


1 


006 





132) 


2.263 





168) 


1.594(0 


111) 


3 


T SCO 


5 


32 





337 





034) 





052 





016) 





728 





060) 





944 





070) 





892 


;o 


075) 





947 





081) 


2.564 





163) 


1.846(0 


107) 


4 


EV Lac 


5 


33 





319 





035) 





038 





013) 





612 





051) 





991 





077) 





845 





072) 





856 





082) 


2.448 





145) 


1.804(0 


101) 


5 


CC Eri 


5 


54 





328 





037) 





038 





014) 





854 





080) 





996 





084) 


1 


011 





088) 


1 


020 





098) 


2.861 





192) 


2.095(0 


126) 


6 


HD 45348 


5 


59 





294 





038) 





039 





013) 





459 





054) 





788 


;o 


076) 





751 


;o 


071) 





959 


;o 


107) 


1.999 





140) 


1.500(0 


108) 


7 


AU Mic 


5 


63 





370 





049) 





025 





014) 





628 





077) 





801 





074) 


1 


008 





096) 


1 


266 





142) 


2.437 





183) 


1.748(0 


124) 


8 


44 Boo 


5 


72 





317 





030) 





038 


;o 


015) 





651 





049) 





916 





064) 





786 





050) 





861 





064) 


2.354 





124) 


1.737(0 


084) 


9 


VW Cep 


5 


72 





409 





047) 





056 





020) 





649 





067) 





927 


;o 


103) 





930 





099) 


1 


013 





135) 


2.507 





192) 


1.711(0 


115) 


10 


Capella 


5 


85 





364 





004) 





056 





002) 





701 





005) 





963 





005) 





896 





005) 





930 





000) 


2.560 





015) 


1.802(0 


009) 


11 


HD 93497 


5 


98 





292 





026) 





086 





015) 





677 





051) 


1 


002 




066) 





872 





067) 





872 





066) 


2.551 





147) 


1.850(0 


093) 


12 


AB Dor 


5 


99 





405 





040) 





059 


;o 


018) 





679 





049) 





927 


;o 


067) 





912 





066) 





987 


;o 


078) 


2.519 





140) 


1.721X?) 


086) 


13 


TZ Crb 


6 


08 





325 





014) 





057 





007) 





603 





021) 





875 





029) 





756 





026) 





866 





035) 


2.233 


;o 


056) 


1.616(0 


036) 


14 


V824 Ara 


6 


17 





332 





030) 





056 


;o 


015) 





615 





047) 





948 





053) 





910 


;o 


054) 





962 





063) 


2.473 





114) 


1.782(0 


079) 


15 


Stack 


6 


19 





345 


;o 


030) 





026 





012) 





640 





044) 





871 





055) 





862 


;o 


050) 





991 





057) 


2.374 


;o 


121) 


1.732(0 


083) 


16 


ER Vul 


6 


28 





344 


;o 


034) 





094 





017) 





598 





053) 





899 





077) 





878 





070) 





980 





080) 


2.375 





161) 


1.651(0 


096) 


17 


HR 1099 


6 


35 





305 





020) 





054 





011) 





626 





034) 





853 





044) 





864 





042) 


1 


016 





059) 


2.343 





088) 


1.725(0 


059) 


18 


/3 Ceti 


6 


42 





347 





016) 





066 





007) 





686 





029) 





929 





032) 





862 





030) 





929 





040) 


2.477 





067) 


1.753(0 


043) 


19 


Algol 


6 


57 





372 





035) 





061 





017) 





686 





049) 





993 


;o 


067) 





857 


;o 


057) 





865 


;o 


065) 


2.536 


;o 


129) 


1.771(0 


083) 


20 


AR Lac 


6 


57 





400 





042) 





014 





011) 





640 


;o 


053) 


1 


017 





070) 





907 




066) 





895 





075) 


2.564 





147) 


1.814(0 


098) 


21 


UX Ari 


6 


63 





314 





050) 





027 





019) 





601 





080) 





631 





088) 





740 





091) 


1 


187 





182) 


1.971 





196) 


1.471(0 


137) 


22 


IM Peg 


6 


88 





327 


;o 


043) 





085 





024) 





610 





063) 





766 





078) 





832 


;o 


077) 


1 


094 





121) 


2.209 


;o 


168) 


1.565(0 


109) 


23 


a Gem 


6 


95 





322 





021) 





053 





012) 





603 





032) 





851 





047) 





883 





047) 


1 


040 





065) 


2.337 





092) 


1.700(0 


069) 


24 


A And 


7 


02 





381 





037) 





059 





019) 





591 





046) 





868 





053) 





793 





057) 





916 





073) 


2.252 





111) 


1.564(0 


073) 
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Fig. 1 . — Comparison of Fe XVII 
( 120031 ); broken lines are those of 



ine ratios in various theories. The solid line are those of Gu 



filled stars are those of 



Loch et al.l (120061 ). open circles are from APEC, fi ll ed tr i 



angles are those of iDoron fe Beharl (120021) . open squares are from lBhatia fe Doschek 



Chen et all J2003h for 3D/3C and 3F/3C ratio, and lchen fe Pradhan 
( 120051 ) for 3s /3d ratio, the open stars are t hose of IChenl (120071) f or the 3D/3C rati o. For 



3F/3C and 3s /3d, the open stars are those of lChen et al.l (20031) and IChen fe Pradhan! (120051 ) 
scaled for the updated 3C and 3D cross sections of IChenl (120071 ). 
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Fig. 2. — Reconstructed DEM of the combined Capella observation. D in the unit is the 
distance of the source in cm. The broken line show the total Fe XVII emissivity weighted 
DEM. 
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Fig. 3. — Chandra HETG spectrum of the combined observation of Capella in the 12-18 A. 
The lower curve is the sum of ±1 orders of MEG data, the upper curve is the model fit with 
the reconstructed DEM and elemental abundances, which is shifted to the upper-right for 
clarity. The Fe XVII lines and the Fe XVI line at 15.21 A are marked. 
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Fig. 4. — Comparison of Fe XVII line ratios calculated with the reconstructed DEM of the 
observations (solid lines), and those calculated at a single temperature of T em (open circles). 
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Fig. 5. — Comparison of measured and predicted Fe X VII line rat ios. Filled circles with error 
bars are the mea sured values, solid lines are those of iGul (120031 ). broke n lines are those o f 



Loch et al.l (120061 ). open circles a re from APEC, , f i lled s tars are those of IChen et al.l (120031 ) 



for 3D/ 3C and 3F/3 C ratio, and lChen fe Pradhanl (120051 ) for 3s/3<i ratio, the open stars are 



those of Chen 



of Chen et al 



2007) for the 3D/3C ratio. For 3F /3C and 3s/3d, the open stars are those 



(20031) and lChen fe Pradhanl (120051 ) scaled for the updated 3C and 3D cross 



sections of IChenl ( 120071 ). 
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Fig. 6. — Comparison of measured Fe XVII line ratios with the pre d iction s of the present 
MBPT corrected theory (solid lines). The predictions of lLoch et al.l (120061 ) are also shown 
as broken curves. 
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Fig. 7. — Comparison of measured and predicted 3C and 3D cross sectio ns. The open circles 
with error bars are the laboratory measurements of iBrown et al.l (12006I ). solid lines wit h no 
symbo ls are the present MBPT corrected calculat ion, dashed lines are the values from iGu 
(120031) . dotted lines are those of lLoch et al.l (120061 ). solid lines with fi l led tr iangles are those 
of I Chen et al.l (120031 ) , and open squares are the calculations of IChenl ( 120071 ) . 
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Fig. 8. — Density dependence of the Fe XVII line emissivity at a temperature of 10 6 7 K. 
The solid l i nes ar e the present MBPT corrected calculations, and dashed lines are those of 
Loch et all tood ). 
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Fig. 9. — Density dependence of the M2/3G line ratio at two marked electron temperatures. 
The solid l i nes ar e the present MBPT corrected calculations, and dashed lines are those of 
Loch et all fhoOfil ). 



